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Asymptotics for a determinant with a confluent hypergeometric 

kernel 



P. Deift^, I. Krasovsky^, and J. Vasilevska^ 



Abstract. We obtain "large gap" asymptotics for a Fredholm determinant with a confluent hy- 
^S ', pergeometric kernel. We also obtain asymptotics for determinants with two types of Bessel kernels 

i^^_^ ' which appeared in random matrix theory. 

(N 

'Th^' 1 Introduction 

r~| , Let K^'^'^' be the operator acting on L^(— s,s), s > 0, with kernel 

-)— » 

^' k(c.,0)(^..^ 1 r{l + a + (3)T{l + a-(])A{u)B{v)-A{v)B{u) 

^ ' ' 2m r(l + 2a)2 u-v ' ^ ' 

cn ■ where 

^ ; A{x) = gl^'^{x)\2x\°'e-'''(p{l + a + /?, 1 + 2a, 2ix), 

^ '. B{x) = 5y^(x)|2x|"e^^0(l + a - /3, 1 + 2a, -2ix), 

O: 5/3(x) = { ^~r/' !^n' > a,f3eC, Ka>-l/2, a ± /3 / -1, -2, . . . 

Here r(x) is Euler's F-function, and (j){a,c,z) is the confluent hypergeometric function (see, e.g., 

^2: </>(a,c,z) = l + y "l" + \^---;" + ^"\^ ^. (2) 

c3 ■ v^v , > ; ^Z^ c c + 1 ••• c + n-1) "I ^^ 

. . . 1 n=l 



n! 



Using the standard recurrence formulae for <j)(a, c, z) (see (36) below), we can rewrite (1) in another 
form:^ 

^ ^ vr (l + 2a)r(l + 2a)2 y/3 v ;y/3 v ; ^_^ ^^^ 

X [n</>(l + a + /?, 2 + 2a, 2iu)(j){a + /3, 2a, 2z7;) - 7;0(1 + a + /3, 2 + 2a, 2iv)(t>{a + /3, 2a, 2in)]. 

The kernel (1) or (3) is called the confluent hypergeometric kernel. For a G M, /3 G iM (in this 
case the kernel is real, which is easy to see from (1)), it was considered by Borodin and Olshanski 
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in [5], and by Borodin and Deift [4] (Proposition 8.13). This kernel arises in several different, but 
related, contexts: 

First, following [5], consider the space H of infinite Hermitian matrices (-f/^jfc)?fc=i- The U{oo), 
the inductive limit of the unitary groups U{N), N — )■ oo, acts on H by conjugations. A probability 
Borel measure on H which is invariant under the action of C/(oo) is called ergodic, if any invariant 
mod set has measure or 1. Consider the space Q whose elements consist of 2 infinite sequences. 



oo 



af > a^ > ■ ■ ■ > 0, a^ >a2 > ■■■ >0, where X](«j')^ + X]("j )^ < °°' 

together with 2 extra real parameters 71, 72, where 72 > 0. It turns out that the elements of il 
parametrize the ergodic measures on H. Furthermore, it can be proved that any C/(oo)-invariant 
probability measure on H decomposes on ergodic components, i.e., it can be written as a continuous 
convex combination of ergodic measures. This spectral decomposition is unique and is determined 
by a probability measure on Q which is called the spectral measure of the original invariant measure. 

The space Q maps to the space Conf(M*) of point configurations on the punctured real line 
M* = M \ {0} in the following way: 

({a+}j=i,2,..., {a7}i=i,2,..., 71,72) -^ (-a^;, -02", . . . , a^,a^), 

where possible zeros among a • are omitted. Under this map spectral measures corresponding to 
invariant measures on H, push- forwards to measures on Conf(M*), give rise in this way to random 
particle systems on M*. 

In [5], the authors considered a particular class of f^(oo)-invariant measures on H, the Hua- 
Pickrell measures, and showed that the push- forwards of these measures to Conf(]R*) give rise 
to random particle systems on R* which are determinantal with correlation kernels given by 

K^°''°'{l/u, l/v)/{uv) (see (1)) with parameters a, i/3 G M. 

The kernel (1) is also (see [4]) a particular scaling limit of a kernel K'^'^''^{u,v) which has a 
similar structure but with the confluent hypergeometric functions replaced by the hypergeometric 
functions 2^1 (ck, /?, 7; z). The kernel K°-'^''^{u, v) is the correlation kernel for a particle system that 
arises in the theory of representations of [/(oo): role of the ergodic measures is now played by the 
indecomposable characters of U{oo) which are again parametrized by certain sequences together 
with some extra parameters (see (1.4) in [4]). 

The kernel (1) also arises as the correlation kernel for a particle system in a similar way to 
i^"'^''^(n, v) above, but in place of irreducible representations of U{oo), we now consider irreducible 
(spherical) representations of U{oo) x H{oo) = lim^ G{N), the inductive limit of the semidirect 
product G{N) = U{N) k H{N), where U{N) is the group of N x N unitary matrices and H{N) 
denotes N x N Hermitian matrices: see [5] and references therein. 

As we will see in Section 2, the kernel (1) can be obtained as a scaling limit in unitary random 
matrix ensembles generated by the weight function f{z, 0) on the unit circle (given by (20) below) 
at a point of so-called Fisher-Hartwig singularity. This singularity combines a root-type and a 
jump-type singularity characterized by the parameters a and /3, respectively. 

In particular cases, the kernel (1) reduces to a Bessel- and the sine- kernel which attracted 
much attention mostly because of their interest for random matrices. If /3 = the confluent 



hypergeometric function reduces to Bessel functions (see, e.g., [1]): 



,^(^,2^,2zx) = r ( ^ + i] e- (f )"''^' J^_i (x). (4) 



Therefore, we obtain from (3) 



\ci\„,\a 



This kernel appeared in [2], [22], [25], [29]. 

If a = 0, /3 = 0, then (5) reduces to the sine kernel 

K('''\x,y)^KUx,y) = "-^^^^^, (6) 

7r(x - y) 

the most ubiquitous object of random matrix theory. 

Note that the operator K^'^'"' is trace class (see Appendix), and consider the Fredholm deter- 
minant 

det(/-K("''3))^2(_,,,). (7) 

Because of the mentioned interpretation of K^'^'^' with a, i/3 G M as the correlation kernel for a 
particle system produced by a Hua-Pickrell measure, it is easy to see that the Fredholm determinant 
(7) is the probability that all the a- are less than 1/s. 

By a random matrix interpretation of the kernel (1), the determinant (7) with a, i/3 G M gives 
the probability, in the bulk scaling limit, that the interval (— s, s) with a Fisher-Hartwig singularity 
at the center contains no eigenvalues of corresponding unitary random matrix ensembles. 

As noticed in [29, 4], the determinant det(/ — -f^^"''^'')L2(o,s) is related to a solution to the 
Painleve V equation.^ 

In this paper, we obtain the asymptotics of the Fredholm determinant (7) for large s, i.e., the 
large gap asymptotics. Our main result is the following. 

Theorem 1 Let K^°''^' be the operator with kernel (1) acting on L?'{—s,s), Then, as s ^- +oo, 



\ )L { s,s) 22"^G(l + a + /3)G(l + a-/3) [ ysj 



, (8) 



where G{x) is Barnes' G-function. This expansion is uniform in compact subsets of the a-half-plane 
Ka > —1/2 and of the [5-plane outside neighborhoods of the points q it /3 = —1, —2, .... 

Remark 2 Setting /3 = in (8) and using a doubling formula for the G-function, we obtain the 
large s-asymptotics for the determinant with kernel (5) where Ka > —1/2: 



det(/ - 41,,,,) = ^G(« + l/2)G(a + 3/2)s-i-^--+2a. 
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1 + 0-1 

s 



(9) 



*More generally, an 2J'i-kernel determinant is expressed [4] in terms of a solution to the Painleve VI equation. 
For some asymptotic results which use this connection and conjectures see [23]. 



Setting a = /? = in (8) and using the property 21nG(l/2) = (1/12) In 2 - ln^+ 3C'(-1), 
where C{x) is Riemann's zeta-function, we reproduce the result for the sine-kernel determinant: 

lndet(/-i^sm) = -y- Jlns + ^ln2 + 3C'(-l) + oQy s ^ oo. (10) 

The first two terms in (10) were first found by des Cloizeaux and Mehta [6], and the full expansion by 
Dyson [14]. The calculations in [6, 14] were not fully rigorous. A proof for the first leading term was 
carried out by Widom [27]. The full asymptotics of the logarithmic derivative (d/ds) lndet(/ — Ksin) 
were proved by Deift, Its and Zhou in [9]. Finally, the constant term in (10) was proved in [15], 
[20], [11]. 

Remark 3 In the present paper, we address only the symmetric case of L^{a,b) such that b = 
—a = s > 0. However, one can apply our methods to consider non-symmetric cases as well. 

In unitary random matrix ensembles at a hard edge of the spectrum (e.g., Jacobi at the edges 
or Laguerre at zero) local correlations between eigenvalues are expressed in terms of the following 
Bessel kernel first considered by Forrester [19] (in an equivalent form): 

^Bessel2{^^y) - 2{x - y) ' ^ ^ 

In particular, the distribution of the extreme eigenvalue is given in the scaling limit by the 

Fredholm determinant det{I — K^^^^^^^) i where -f^)^essei2 ^^ ^^^ trace-class operator on L^(0, s), s > 0, 
with kernel (11). In Section 7, we prove the following asymptotic behavior of this determinant. 

Theorem 4. As s ^ +oo, we have uniformly in compact subsets of the half-plane Ka > —1.' 

det(/ - 4tse^2)L^(o,s) = r,.-«V4e-/4+aVi ^^ + 0{s^'/')) , Ka > -1, (12) 

where 



G{1 + g) 

(27r)"/2 



(13) 



.(a) 



In [26], Tracy and Widom showed that the logarithmic derivative (d/ds) In det{I — K^^^^^^^) 
is expressed in terms of a solution to Painleve V equation and used this fact to give a heuristic 
derivation of (12) with some constant Ta- Tracy and Widom also conjectured the value of Ta given 
in (13) using numerical calculations and comparison with the Dyson asymptotics for the sine- kernel 
determinants. (In fact, for a = =Fl/2, the Bessel kernel (11) reduces to sine-kernels appearing in 
orthogonal and symplectic ensembles of random matrices. The sine-kernel (6) appears in unitary 
ensembles.) Very recently, a proof of the asymptotics (12,13) for the range of the parameter |Ka| < 1 
was given by Ehrhardt [16] using operator theory methods. 

To prove Theorem 1, we use the approach of [20], [11], [12], where the asymptotics were 
computed, including the constant terms, of the sine-kernel and the Airy-kernel determinants. 

First, in Section 2, using results from [10] we express (see Lemma 6) the Fredholm determinant 
(7) as a scaling limit of Toeplitz determinants Z)„(v) with certain symbols /(e*^) supported on an 



arc of the unit circle (p < 6 < 2it — (p with p = 2s/n, n > s. The continuation of these symbols 
into the complex plane has a Fisher-Hartwig singularity at z = 1. Theorem 1 then reduces to an 
asymptotic evaluation of such Toeplitz determinants for large s. 

In Section 3 we derive a differential identity (49) for the logarithmic derivative {(P' /dip^) In Dn{p) 
at < c/? < vr in terms of the solution to an associated Riemann-Hilbert problem (in fact, in terms 
of the associated orthogonal polynomials which are given by this solution). In section 4, we obtain 
the series expansion of Dn{<p) for (p close to tt. In Section 5, we solve the Riemann-Hilbert problem 
asymptotically and thus obtain the asymptotic expression for the r.h.s. of the differential identity 
(49), namely, we obtain the identity (139) uniformly for 2s/n < p < n, n > s, s > sq, with some 
(large) sq > 0. Integration of the latter identity w.r.t. (p, using the boundary condition of Section 
4, gives the asymptotics of the determinants Dn{<p) for any arc with 2s/n < p < n, n > s, s > sq, 
with some sq > 0, which is sufficient to prove Theorem 1. 

In Section 7, we represent the Fredholm determinant with the Bessel kernel (11) as a scaling 
limit of Hankel determinants related, via a general connection formula of Theorem 2.6. of [10], to 
the particular case of Dn{p) with /5 = 0. The connection formula also involves the polynomials 
orthogonal w.r.t. f{z) on the circular arc which are represented by matrix elements of the solution 
to the Riemann-Hilbert problem mentioned above. We prove Theorem 4 by using asymptotic results 
on these polynomials and on Dn{p) from the previous section as well as an expansion for singular 
Hankel determinants from [10]. 

2 Connection with Toeplitz determinants 

The aim of this section is to derive an expression for (7) in terms of Toeplitz determinants (Lemma 
6 below) and to fix notation for the rest of the paper. 

Let E^ be an arc of the unit circle C oriented counterclockwise: 

E^ = {e'\ p<e<2T:-p], < 99 < vr. (14) 

Consider the following function /(z, p) on Ei^: 

f{z,p) = \z-l\^'^z^e-'^l^, z = e'^eE^ a,f3eC, Ka > -^. (15) 

Note that for z G C 

l^_i|2«=(i^l)!!, (16) 

where the cut of {z — 1)^" is along [l,oo), and < arg(2; — 1) < 2tt. The branches of z", z^ are 
chosen so that < argz < 27r. Therefore, we can extend the function f{z) to the complex plane 
with the cut [0, 00) by the expression: 

/(^) = z-"+'^(z - l)2°e-^^("+^) zGC\[0,oo). (17) 

Related to the function (15) is a system of orthogonal polynomials pk{z]Lp) = Xki^)^^ + . . ., 
Pk{z; ip) = Xk{^)z'' + ... of degree A: = 0, 1, . . ., satisfying 

27r— (p 2n—ip 

^ J Pk{z)z-^f{z)d0 = x;n'Sknu ^ J Pkiz~')z"'f{z)de = X-^r^hm, (18) 



z = e , m = 0, 1,. . . ,/c. 

Note that if the weight function f{z) is not positive on E'^, the existence of such a system of 
polynomials is not a priori clear and will be addressed in the situations needed below. 

In order to obtain the kernel (1) in a scaling limit, we will need to know the asymptotics of 
the polynomials 

g„(z) =p„(z;0), g„(z) =p„(z;0), (19) 

corresponding to the weight 

/(z, 0) = |z - l|2°z/^e-*^^, z = e\ < e < 27r. (20) 

In this case the function /(z) possesses a Fisher-Hartwig singularity at the point z = \. The 
asymptotics of these and more general polynomials were recently analyzed in [10] (in particular, 
the polynomials exist for sufficiently large degrees) and from those results we obtain the behavior 
of the polynomials in a neighborhood of the singular point. 

Lemma 5 Let < e < 1, f^o = {2;, |z — 1| < e\. Fix the branch o/lnz = ln|z| +iargz hy the 
condition — vr < argz < vr, and the branches of the power functions w'^ = |w|"exp{zaarg w} by the 
condition < argw < 27r. Then, as n ^ 00, z G Uq, 

/ N n (21) 

x^^7^-^-^0(l + a + /3,l + 2a,nlnz) l + o(- 
1 (1 + 2a) I \n 

g„(z-i) = (nlnz)"+'^(z - i)-("+/3) ^(^ + "~(^) ^(i +a-P,l + 2a, -nlnz) 

r(l + 2a) 



n 



(22) 



These asymptotics are uniform and dijjerentiable for z £ Uq. They are also uniform in compact 
subsets of the a-half-plane "^a > —1/2 and of the f3-plane outside neighborhoods of the points 
a±p = -1,-2,.... 

Remark It is easy to check that the singularities in (21), (22) cancel. 

Proof. Our polynomials correspond to the special case of [10] with only one singularity located at 
z = 1. Let Uq be the neighborhood of 1 where the parametrix was constructed in [10] in terms of 
the confluent hyper geometric function. Let 

Q = nlnz. 

Take ( £ I, where / is the first sector of the image of the neighborhood under the conformal 
transformation ^ = nlnz (see Figure 2 of [10]). Tracing back the Riemann-Hilbert transformations 
of [10], it is straightforward to obtain 






n-1^ , ,, ,-,/ + n--/^-0 5 n-/^- -M, n^oo (23) 



where 



r(l + a + 0) (24) 

T{a - /3) 
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*2 = ^l^{a + |3^ + 2a, () - ^(1 + a - /3, 1 + 2a, e— C)^^^^^^^^^-^-"^". ^^^^ 

i (a + p j 

Here 'i}j{a,c,z) is the confluent hypergeometric function of the second kind (see, e.g., [1]), and 
0{l/n) stands for a 2 x 2 matrix with the matrix elements of that order. 

Applying the following property of the confluent hypergeometric functions: 

-p/-i \ -p/ -|\ 

Ha, c, z) = '"^^. Ha, c, z) + l^" ' z^-'<p{a -c + l,2-c,z). (26) 

i (a — c + Ij I (a) 

and Kummer's transformation 

4i{a,c, z) = e^<l){c — a,c,—z) (27) 

to ^2 in (25) gives 

^' = ^r(l^+2a)^^ ""^^^"^^^^"'^^^ + a - /3, 1 + 2a, -C), (28) 

which simplifies the expression for Q2. To simplify the formula for Qi, we use (26) and (27) again. 
We obtain for the combination ^i in Qi: 

^1 = e'-^^^^f^e-'-^-^^Uil + a + /3, 1 + 2a, C). 

Now the expressions (21), (22) of the lemma for C G / follow easily (noting also that Xn = l+0{l/n) 
by Theorem 1.8 in [10]). As Quiz), Qn{z) are polynomials and the asymptotics of [10] hold uniformly 
for z G C/q) the expressions for quiz), Qniz) extend by continuity to the whole neighborhood Uq 
and hold there uniformly. The uniformity properties in a and /3 follow from the uniformity of the 
asymptotics in [10]. The multiplier e-27rj(Q!-/3) ^^^ qn{z) in C_ n Uq appears because of the cut of 
z"- going through the neighborhood. D 

Let Dn{^) denote the Toeplitz determinant with symbol f{z,ip): 

Dn{^)=dei{f,_,)-~l, = -^l' ...j n \z,-z,\'flf{z„^)^, (29) 



E^ E^ i<i<fe<« i=i ■> 



where ft are the Fourier coefficients of /(z, ip): 

-I /'27r— tp 



fk = ^J^ ^ fi^'', V)e^"''d0, k = 0, ±1, ±2, . . . 



Then the following lemma holds. (For a = /3 = it reduces to the scaling limit used by Dyson 
in his analysis of the sine-kernel determinant [14].) 

Lemma 6. Let s > 0. Then 



D„ m 



det(/-i^(°''^))= lim " \"/ , (30) 

n^oo D„(0) 



where K^^'"' is the operator on Lp'{—s,s) with kernel (1). 

Proof. Assume first that a,i/3 € M. Then, as follows, e.g., from (29), Dn{(p) > for all n, 
and therefore the polynomials qk{z), qk{z) exist for all k (as follows from their determinantal 

representation: see, e.g., [7]). By a standard argument [24, 7], we first write the term W \zj — z^p 

j<k 
in (29) as a product of two Vandermonde determinants whose elements, by a suitable combination 
of the rows, become the polynomials (19) qk-i{zj)/xk-i{0) and qk-i{z~ )/xfc-i(0), j,k = 1, . . . ,n, 
respectively. We obtain 



j<k 



n-l 



/n— 1 



" - n ^^(0)"' ^ det lY.^,izj')qeizk) 

j=0 \e=0 / l<j,k<n 

Using the well-known expression 



we obtain 



1 



Dn{0) ~ (27r)"n 



n-l 
j=0 



det{Kn{zi,Zj))^^.^.^^ 



dzi dZn 



IZl 



where the kernel Kn is given by the expression: 






(31) 



rt— 1 



Kn{zi,Z2) = ^/f{zi,0)f{z2,0) ^ gfc(^i ^)qk{z2) 



fc=0 



Vf{zi,0)f{z2,0) 



[z2/ziYqn{zi)qn{z2 ^) - gn(^i ^)gn(^2) 
1 - Z2/Z1 



(32) 



To obtain the second equality here, we used the Christoffel-Darboux formula: see, e.g.. Lemma 2.3. 
in [10]. Since for sufficiently large n, both the polynomials quiz), qn{z) exist (see Lemma 5) and 
-^n(O) / (see (38) below) for complex a and /3, equation (31) is extended to the general case of 
a, (3 from a, i/3 G M by continuity and holds for all sufficiently large n for a and /3 in a compact 
set. 

As, e.g., in [7] Section 5.4., one shows that the r.h.s. of (31) can be written as the Fredholm 
determinant det(/ — Kn), where Kn is the operator with kernel (32) acting on L^(C \ £^^, 2^)'-' 
that is 

Dni'p) = Z?n(0) det(/ - Kn)L^i^c\E^), (33) 

where the arc C \ E^ is oriented counterclockwise. 

We now show that the kernel (1) can be obtained as a scaling limit of (32). Setting z = e~n~ , 
u > 0, in (21), (22), we obtain 



g„(e^) = ^""^ ^r(]^^2t)^^ 0(1 + a + /3, 1 + 2a, 2ni) 
Ue-'^)=n-+P^±^^<P{l + a- P^ + 2a,-2^u) 



i + o[- 

n 



1 + 



n 



(34) 
(35) 



Setting now z = e '^^^ « , u < 0, in (21), (22), we obtain the same expressions as the r.h.s. of (34), 
(35) for the values of the polynomials q^ Qn at these points. 

Let z{w) = e~^ if it; > 0, and z{w) = e ^^~^~^ if w < 0. Substituting the just found values of 
the polynomials into (32), we obtain 



lim Kn{z{u),z{v))^^^^ = e^(^-")i^("'^)(n,f)(it;, u,v e 



2TTiz{v) ' ' ' ' 2TTiz{v) vrn' 

where K^°-'^'{u,v) is given by (1). 

Using the following standard recurrence relations for the confluent hypergeometric function 

a(j){a + 1, c, x) — (a — c + l)4>{a, c, x) — (c — l)0(a, c — 1, x) = 0, 
c4){a + 1, c, x) — c(t){a, c, x) — x(j){a + 1, c + 1, x) = 0, 

we can rewrite (1) in the form (3). 

Now using the estimates (34)-(35) for the polynomials qniz), we see that for any s > there 
exists c(s) > such that 



did^ (—Kn {zi,Z2) - e^(-"")i^("'^)(u,«)) 
\7rn / 



< -, (37) 

n 



where u,v £ (— s, s), j, k = 0,1. By similar arguments to those of the proof of Corollary 1.3 in [8], 
the estimate (37) leads to (30). D 

A well-known result on asymptotics of Toeplitz determinants with a single Fisher-Hartwig 
singularity such that a ± /3 / —1, —2, . . . (see, e.g., [17, 10]) reads in our case: 

Z,.(„) = „<.'-^^^l(l±i!±»ii^(l+„(l),. (38) 

This expansion is uniform [10] in compact subsets of the a-half-plane Ka > —1/2 and of the /3-plane 
outside neighborhoods of the points a it /3 = —1, —2, .... In the next sections we will obtain an 
expression for Dn(2s/n) for large s, n, s < n. 

3 Riemann-Hilbert problem and a differential identity 

Suppose that the system of orthonormal polynomials satisfying (18) exists and consider the following 
matrix-valued function: 

/ -1 c^^ -1 r PnMimdL \ 

^^"ly iz«-ip i(z-i) y if 2^.1^1im^]' ^^^^' ^^^^ 

where C^, is the arc E^p (14) but oriented clockwise. Denote z+ = e**^, Z- = e^^'^^~^' , the endpoints 
of the arc. 

It is easy to verify directly that Y{z) = y'")(^) solves the following Riemann-Hilbert problem: 



(a) Y{z) is analytic for z £ C\ C^. 

(b) Y{z) has continuous boundary values l + (-z) as z approaches the inner points of the arc C^ 
from the outside of the unit circle, and Y^{z), from the inside. They are related by the jump 
condition 

y+(z) = y_(^)(j ^^Y^^^)' ^eG^\W,^-}- (40) 

(c) Y{z) has the following asymptotic behavior as ^ — )■ oo: 

Y{z)=[l + o{^\\z^''^\ where as = (J _M . (41) 

(d) Near the endpoints of the arc, 

n'^-o{\ I:!::::!). <-> 

as z — > z-t, 2: G C \ C^p. 

The solution (39) to the RHP (a)-(d) is unique. Note first that dety(z) = 1. Indeed, from 
the conditions on Y{z), det y(2;) is analytic across C^^ has all singularities removable, and tends to 
1 as 2: — >• 00. It is then identically 1 by Liouville's theorem. Now if there is another solution Y{z)^ 
we easily obtain by Liouville's theorem that Y {z)Y {z)^"^ = 1. 

A general fact that orthogonal polynomials can be so represented as a solution of a Riemann- 
Hilbert problem was noticed for polynomials on the real line by Fokas, Its, Kitaev in [18], and 
extended to polynomials on the circle in [3]. The point of this representation is that the Riemann- 
Hilbert problem can be efficiently analyzed for large n by a steepest-descent method discovered 
by Deift and Zhou [13] (and developed further in many subsequent works). This gives the large-n 
asymptotics of Y(z), and therefore, by (39), the asymptotics of the orthogonal polynomials. We 
defer the asymptotic analysis in the present case to Section 5. The rest of this section will be 
devoted to a derivation of a differential identity for D„(c/?) in terms of the matrix elements of Y{z). 

We start with the following auxiliary lemma (which, in fact, is true for any weight f{z) and a 
jump contour C^). 

Lemma 7. Let the system of polynomials Pk{z), Pk{z), k = 0, . . . satisfying (18) exist. Fix n > 1. 
Then we have the following Christoffel-Darhoux identity in terms of the function (39): 



^M^'')pu{z) = -z-"+i limtr (^ (° J) y(C)-i) , UC^- (43) 






Remark. The right hand side of (43) contains only the elements of the first column of Y which 
are analytic. 

Proof. Multiplying the recurrence relation (2.4) of Lemma 2.1. in [10] by z, and replacing n with 
n — \ gives (in the present notation) 

ZPn[.Z ) = Pn-l{z ) H Z ^ Pn{z)- 

Xn Xn 
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Substituting this expression into the r.h.s. of the Christoffel-Darboux identity (2.8) of [10]: 

n-l 



^Pk{z ^)pk{z) = -npn{z)pn{z ^) + z(p„(z ^)-T-Pn{z) - Pn{z)—Pn{z ^) j , 



k=Q 



-ud 



d 



we obtain 



n— 1 



EPk{z ^)pk{z) = -{n-l)— Z ^Pn{z)pn-l{z ^) 

fc=0 ^" 



Xn-l 
Xn 



Pn{z)—pn-l{z ^)-pn-l{z ^)-j-Pn{z) 



(44) 



(n - l)z-"y2i(^) + z~'^+'^£Y2i{z), we obtain from (44) 



Using the expressions Yuiz) = Xn^Pn{z), Y21 = Xn-i^" ^Pn-i{z ^), and Xn~ij^Pn~i{z ^) 

A. 

dz 

n-l 

n,Jz-^)n,Jz) = z~''+'^(Yo, .,, .,, 

dz dz 



Y.Pk{z-^)Pk{z) = z~'^+\Y2i^Yn - Yn^Y2i), 

fe=0 



(45) 



which proves the Lemma. □. 

We are now ready to formulate the main result of this section. 

Lemma 8. Let the polynomials Pk{z), Pk{z), k = No,No + 1,..., satisfying (18) exist for some 
Nq > 0. Fix n > Nq. Let z^, z_ be the endpoints of the arc C^: 2+ = e**^, z_ = e*'-^'^^'^-*. Then 



li_lnZ)„(v9) = trCo(C+ - C_) + ^+-trCiC+ - 

I dip zj^ — 1 Z- — \ 



^" trCiC_ + ^+^^~ trC+C_, (46) 

z I — z_ 



where, in terms of the matrix (39), 

n + a - (3 
C/n = ; 



y(0)^3^"'(0), Ci = aY{l)a3Y-\l), 



(47) 



and the following limits taken for z ^ dp exist and define C± : 



^-=.'isv^^'w(2 >-'^'^- ^- 



A-?_^^'w(s ;i^'-'w. (48) 



with the extension of f given by (17). Moreover, the second derivative 



dip' 



■InDni'p) 



^+ trCi(C+ + C-) + — -tiC+C-. 



{z^ - 1)2 



U+ - -z- 



(49) 



Proof. Assume first that a, i/3 G M. Then we can set A'g = (cf. proof of Lemma 6). Starting with 
the representation of a Toeplitz determinant in terms of the leading coefficients of the polynomials 
Pk'- 



n-l 



Bn{^) = WXk\^). 



(50) 



fc=0 
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we obtain, using first (18) and then integration by parts, that 



J ' { \ 1 pZlT — 'p p) 



^^X.(V^) 2W^ f^, 

^ n— 1 



(51) 



A;=0 

The previous lemma immediately gives 



d ^ ^ ^ ^ -n+l,. .. r.^ .JdY..(Q Mw.^-l 



2vr-lnl)„(^) = .;-V(.„^)^Hm tr ^-(.) ^^ ^ ) y(z) 

+._-v(.-,..)^Hmtr(f (.)(;; ;)y(.)-). 



(52) 



In the next section we will show that the solution Y of the Riemann-Hilbert problem exists for all 
n > Nq with A'^o sufficiently large uniformly for a and /? in a compact set. Therefore, the identity 
(52) extends to the general complex a and (3 from a, i/3 G M by continuity. 

One could already use the identity (52) for the purposes of the present paper. However, 
following the philosophy of the Riemann-Hilbert-problem approach [9], we can simplify it further 
to the form (46) which does not contain derivatives of Y . In order to do this, consider the function 

Y{z) = YizMz^/^ uiz) = z'^f{z), (53) 

where / outside the arc is given by (17). The function Y{z) is easily seen to be the solution of the 
problem: 

(a) Y{z) is analytic for z G C \ (C^ U R+); 

(b) On the contours C^p and M^^: 

y+(z) = y_(z)(j I), zGC^\{z+,z_}; (54) 



y+(z)=y_(z)(^) , zG(0,+oo); (55) 



(c) Y{z) = {I + 0{\)) z'^fizf-i, z ^ oo. 
Since —^ is constant on (0,+oo), we see that 

dY 



^<^> - dT^" 



has no jumps. Since 
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we obtain from the condition (c) for Y that F{z) = 0{l/z) as z — )■ oo. The Riemann-Hilbert 
problem shows that this function can have isolated singularities at z+, 2_, 0, and 1. First, we 
obtain using (56): 

F{z) = ^+To{z), z^O; F{z) = -^ + T,{z), z ^ 1, (57) 

z z — 1 

where Co, Ci are given by (47), and Tj{z) are Taylor series. 

Now let f7 be a disk of a sufficiently small radius centered at z^. If y is defined by the 
expression 

y(z) = y(z)(J 2^M^-^+)^, ,^u, (58) 

then it follows from (54) that Y{z) in U has no jump, and from (42), that its singularity at z+ is 
removable. Thus Y{z) is analytic in U. Using (58) we then obtain 

F{z) = -^ + Tsiz), z ^ z+; (59) 

Z — Zo- 



1 .^,. . fO I'w-^i^ 



where T^^z) is a Taylor series and 

Using the definitions (58) and (53), we obtain the expression (48) for C+. Note that the limit in 
(48) exists as the logarithmic singularity of Y{z) at z+ cancels from that expression. 

A similar analysis at z_ gives that 

F(z) = -^— +r4(z), z^z_; (61) 

z — z^ 

where T4(z) is again a Taylor series and (7_ is defined in (48). 

Thus we conclude that -F(z) is a meromorphic function with first-order poles at 0, 1, z+, z_, 
and since -F(z) = o(l) at infinity, we have identically in the complex plane 

F{z) = — + — ^ + ^- + 



z z— 1 z — z^ z — z^'' 
or recalling the definitions of -F(z) and ^(z). 



f^(^, V^) = f ^ + ^ + ^ + ^- ) 1^(^, 9^), (62) 

az \z z — 1 z — z+ z — z_ 



Substituting (63) into (52) and noticing that C\. = gives the identity (46). 
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To obtain the identity for the second derivative, note first that as follows from the general 
theory the function Y{z) is differentiable w.r.t. (p. Similarly to our derivation of (62), we obtain 

Equating the derivatives 4- -j-Y {z , (p) = 4^ -j-Y {z , ip) , gives by (62) and (64) a compatibility con- 
dition on Co, Ci, C+, C_, and their derivatives w.r.t. (p. Equating the coefficients at 1/ z in this 
condition gives 

-^Co + i[Co,C+-C_] = 0, (65) 

dip 

where [^, -B] = AB — BA. Similarly the coefficients at l/(z — 1), l/(z — z+), l/(z — z_) yield the 
identities: 






Z^—\ Z_ — 1 



0, (66) 



-f C+ - i[Co,C+] - ^£+-[Ci,C+] +z^i±^[C+,C_] = 0, (67) 

dip Z-^ — 1 Z-\- — Z— 

— C_ +i[Co,C_] + — [Ci,C_] - i^^-[C+,C_] = 0. (68) 

dip Z- — \ ^+ — -2- 

Differentiating (46) w.r.t. ip and substituting the above identities for the derivatives of Co, Ci, 
C+, C_ in the resulting expression gives the formula (49). In this calculation, it is convenient to 
use the elementary algebraic identity: 

tT[A,B]C = iTA[B,C], (69) 

which equals if any two of ^, B, C coincide. D 

4 Expansion of Dn{^) as (/9 — > tt. 

For a fixed n > 1 we will now obtain the expansion of Dn{(p) as y? — t- tt. We use the representation 
(29) of Dn{(p) as a multiple integral: 

Dn{^) = W^r~^ ■ ■ ■ r~^ n \e''^-e''^\'fl\e^'^-l\''^e^'^'e-^-^de,. 

( T^) n. J^ J^ l<j<ck<n j=l 

For the analysis of Dn{p>) as c/? — )■ vr, set (/? = vr — e, e > 0. Substituting 9j = tt + exj in the 
integrals, we obtain: 



An + 0(£^)) , as £ — ^ 0, n fixed. 



(27r)" 
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where 

is a Selberg integral (or a product J][fc=o -^k ' '^^ere x^ are the leading coefficients of the orthonor- 
mal Legendre polynomials). The error term in (70) is of order e^. Indeed, it is easy to see that the 
expansion of the factors with the absolute value in the integrand in (70) gives an error of order e^. 
The factors e^^^i^ produce the following term of order e: 



^ / • • • / 2^ ^i n (^i - ^fc)' n '^^. 



which is equal to 0, as the change of variables Xj — )■ —Xj shows. Therefore, the error term in (70) 
is indeed O(e^). 

The asymptotics of An as n — )■ oo are [28] : 

In An = -n^ln2 + nln(27r) - -InnH In2 + 3C'(-1) + o(l), n -^ oo, 

where Ci^) is the derivative of Riemann's zeta- function. Therefore, 

lnDn{(p) = n2ln(7r-99) + (2Qn-n2)ln2--lnn + — ln2 + 3C'(-l) + (5„ + 0„(e2), e = -ir-ip. (72) 
Here 6n depends only on n and (5„ — t- as n — )■ oo. The term 0„(e^) — ?■ 0, as e — )■ 0, n fixed. 

5 Asymptotic analysis of Y(z) 

We now analyze the Riemann-Hilbert (RH) problem of Section 3 for Y{z) in the limit of large n. 
The analysis is similar to that of [20] and [11]. Consider the function 

^(z) = ^ (z + 1 + y/iz-z+)iz-z-)) , 7 = cos((^/2), (73) 

which conformally maps the outside of the arc C^ onto the outside of the unit circle. Note that 
'I'-l-'I'_ = z for z G C^ \ {z^, z_}. Furthermore, we see as in [20] that 

< 1, for \z\ / 1. (74) 



M/(z) 

We apply several transformations to the Riemann-Hilbert problem. First, set 

T{z) = 7-"'"»y(z)^(z)-"'^^ (75) 

Then we obtain a RH problem which is normalized to / at infinity: 
(a) T{z) is analytic for z G C \ C^. 
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-1 



1 tan 



-itan^- 



Figure 1: Confornial mapping. 



(b) T(z) has L^ boundary values on C^p related by the condition 



for z eC^\ {z+,z_}. 



(76) 



(c) 



T{z) = I + 



as z 



oo. 



(77) 



As in [11], we now go over to the variable A given by the following linear- fractional transforma- 
tion that maps the arc C^p onto the interval [—1, 1] with the point z = z- corresponding to A = — 1, 
and z = Z-! , to A = 1: 



A 



Z + l if 

ztan — , 

z- 1 2' 



\ + i tan ^ 



A 



itan ^ 



(78) 



The complementary to C^p arc of the unit circle is mapped to M \ [—1, 1]. The points z = 0, 1, oo 
are mapped to A = — itan ^, oo, itan ^, respectively. The cut of the function f{z), (0, 1)U (1, -|-oo) 
becomes (— itan ^, —ico) U (ioo, itan ^) in the A-plane (see Figure 1). 

For the case of a varying arc when c/? — )■ and c/? — t- tt, in the RH analysis, one would need to 
consider contracting neighborhoods of the end points z± in the z-plane (cf. [20]). This could be 
carried out. However, in the A-plane, we can keep neighborhoods of the points A = ±1 fixed, which 
considerably simplifies the calculations below. Thus, going over to A is not essential, but useful. 

The problem for T corresponds to the following one in the A-plane: 

(a) T(A) is analytic for A G C \ [—1,1]. 

(b) The boundary values of T(A) on (—1, 1) are related by the condition 



r+(A)=T_(A) 



^+(A) 




-2n 



/(^(A)) 
$_(A) 



-2n 1 



for AG (-1,1), $(A) 



z(A)V2 • 



(79) 
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Figure 2: Contour for Riemann-Hilbert problems. 



(c) 



r(A) = 1 + 



as A 



oo. 



The solution T(A) is related to T{z) by the expression: 

T{z) = Tof{X{z)), 



Tq = T (i tan — ■ 



For the function (17) we have in the A-variable (we denote it /(A) for simplicity): 



/(A) ^ fiziX)) 
The function <&(A): 



-J7r/3 



2 tan 



"P 



2a 



X — i tan 



"P 



-a-l3 



X + i tan 



P 



-a+l3 



$(A) 



X + isin^Vl^^ 



cosf(A2+tan2f)V2- 
Note that, because of the properties of '^{z) discussed above, we have 

|$(A)|>1, A^R; |$±(A)| = 1, AG [-1,1]. 



(80) 



(81) 



(82) 



(83) 



(84) 



Following the steepest-descent method of Deift and Zhou [13] we now change the RH problem 
so that the oscillating behavior of the matrix elements in (79) is converted into the exponential 
decay as n — )• oo. 

Namely, consider the system of contours shown in Figure 2. Let / be the region bounded by 
the curves Si and S2 = (—1, 1); region // is the one bounded by the curves S2 and S3; region III 
is the rest of the complex plane. Define the function S{X) as follows: 
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in region /, 



in region //, 



S{X) - T(A) (^ _j,^^yi^^^y2n ^ j > (85) 

in region III, 

S{X) = f{X). (87) 

The Riemann-Hilbert problem for S is then the fohowing: 
(a,b) S{X) is analytic in C \ (Si U S2 U S3) with the following jump conditions on the contours: 

5+(A) = 5_(A)(^^^^_lJ^^^^_2n J), AGS1US3, (88) 

5+(A) = 5_(A)(_^°^^_i ^{3-^^), AGS2^(-1,1). (89) 

(c) As A — )■ 00, 

^(A) = / + oQ). (90) 

For S to have these properties, the contours Si^3 should not intersect the real axis and the 
cuts (— ioo, —i tan 99/2), {itanip/2,ioo) of /(A). 

Below we will investigate the inequality in (84) in more detail and will show that Si^3 can 
be chosen so that for cp satisfying 2s/n < ip < n, n > s, s > sq, the jump matrix on Si U S3 is 
uniformly close to the identity up to an error of order e~^^^, e > 0, outside neighborhoods of the 
endpoints of the arc. The error is small for sq sufficiently large. This suggests that outside some 
5- neighborhoods Us, Us of the endpoints, the function S can be approximated by a parametrix 
which has a jump only on S2. The problem for this parametrix in the outside-the-neighborhoods 
region is standard [7] and will be presented below. It is solved explicitly. Then we will consider the 
neighborhoods Us, Us, and construct (following [21]) local parametrices there in terms of Bessel 
functions. We then match the outside and the local parametrices on the boundaries dUs, dUg for 
large n, which produces the asymptotic expansion of S in the inverse powers of nsin ^. The latter 
expression is large provided again that 2s/n < (p < n, n > s, s > sq, and sq is sufficiently large. 

5.1 Outside parametrix 

The parametrix outside neighborhoods of A = ±1 is the solution to the following RH problem: 

(a) N{X) is analytic for A G C \ [-1, 1], 

(b) A^(A) has L^ boundary values A^_|_, N_ on (—1, 1) related as follows: 

7V+(A) = iV_(A)(_^°^^_i ^^^^y AG (-1,1), (91) 
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(c) 



N{X) = I + Ol-] , asX^oo. 
As is easy to verify, this problem has the following solution: 



(92) 



N{X) = l{V^r( .l'^^' 

2 \zia — a ) a + a 



-i(a — a ^) 



v{\)- 



a{\) 



A-1 

A + 1 



1/4 



(93) 



where the branch of the root is chosen so that a(A) — )■ 1 as A — t- oo. The Szego function 'D(A) is the 
solution to the following RH conditions: a) T){\) is analytic in C \ [—1, 1]; b) D-|_(A)D_(A) = /(A) 
for A G (—1, 1); c) V{X) — )• const as A — )■ cx). We have 



V{\) = exp 



^/r^A? fi ln/(r/) d'n 



2Tri 



1 v^i-7?2??-A y ' 



(94) 



with the integration over the upper ("+") side of the interval (—1, 1): we choose ^/x > for x > 
and < arg(A ± 1) < 27r. Finally, 



Poo = hm V{X) = exp -- / —= 
A->oo y2TT J_i ^\ 



ln/(r?)dT?\ 
^2 J 



(95) 



In what follows, we will need an expansion of D(X) at the endpoints ±1. The integral in (94) 
can be written as half the integral around a loop encircling [—1, 1]. Deforming the loop, we obtain 



TTl 



^ ln/(?7) dr] „„ -i ./xx , 



7ri lim 



(a + /3) 



itan - 



iR 



idri 1 



idf] 



{^2 _ 1)1/2 n-\ 



(96) 



where — i(x^ — 1)^'^ = Vl — a^^ > on the upper side of (—1, 1). Expanding —^ near A = ±1, we 



obtain: 



■q—\ 



itan - 



iR 



idr] 



1 



(^2 _ 1)1/2 ^_^ 

idf] 1 



-.ij (r?2 - 1)1/2 ^ _ A 
Substituting (97), (98) into (96), we obtain for (94) 



±z-e^^^/2 + 0(A^l), A^±l; 



(97) 
(98) 



P(A) = /(l)V2exp Pi(l) 



A-1 



1/2 



+ ^^2(1) 



A-1 



3/2 



+ 0(A-1 



,5/2 



A^ 1, 



P(A) = /(-1)V2 exp ( Pi(-l) ('^"l '^' + P,(-l) (^Y' + O (A + 1)^/^^ 



2e»- 



2e*' 



(99) 



A 



-1, 
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where 

Pi(±l) = 2(a(l-sm-) ±il3 cos-) , (100) 

and the exact value of T>2i^^) wiU not be used as it cancels from the final expressions below. It is 
clear from the construction that this expansion is uniform in f as well as in q, /? in a compact set. 

5.2 Local parametrices 

Let Us and Us denote the (nonintersecting) 5- neighborhoods of the points 1 and —1, respectively: 
see Figure 2. We choose S to be sufficiently small, see below. 

We will now write down essentially known (see [9, 21, 20]) parametrices P, and P in Us and Us, 
respectively. The parametrices have the same jumps as S inside these neighborhoods and match 
A^ on the boundaries dUs and dUs to the leading order. 

Consider the function 

w(A) = ln2$(A), 

which, for a sufficiently small 5, is analytic inside Us and maps it conformally onto a neighborhood 
of zero. It has the following expansion at 1: 



^(A) = 2nsin2||l-^ 



5 

cos ifi + - 
4 



u + Oiu'^)}, u = X-l, XgUs, (101) 



and 

V^= \/2ni/2sin^(l + 0(n)), n = A - 1. 

This expansion is uniform in 99. 

Consider the following mapping of Us'. 

C = n^uj{X), X£Us. (102) 

For our analysis below, we need \(\ to be uniformly large in f and A on the boundary dUs- 
We see from (101) that this is indeed so if 99 satisfies the condition: — < ip < t:, n > s, s > sq, 
with So sufficiently large. 

The local parametrix in Us is given by the following expression (cf. [9, 21, 20]): 

P(A) = ^(A)Q(n^a;(A))e~"V^^"V(A)""^/^ A G Us, (103) 

where 

E{X) = l=N{X)f{Xr/' ( ^ 7 ) (^^V^)r^', (104) 

and the function Q{C) is expressed in terms of modified Bessel and Hankel functions: 
1) in the intersection of region / in the A plane and Us 

1 / i7i')(e--/2^i/2) H^^\e~'^/\y^) \ 

^^^^ = 2 |^<l/2 (F«)'(e--/2CV2) ^^1/2 ^^(2)y (^^^2^1/2) j ' (105) 
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2) region // and Us 






1 



^(l)(gW2^1/2) 



2 1-^1/2^^2) (gW2^1/2) T,(^l/2lH, 



QiO 



r{l) 



t/2^1/2 



(106) 



(107) 



3) region /// and Us 

,vrzC^/2/^(CV2) _^i/2^/(^i/2);' 
where — vr < arg((") < vr. 

Using the asymptotic expansions of Bessel and Hankel functions for a large argument, one 
obtains uniformly on dUs'- 



P{\)N{\)-^ =1 + N{X)f{Xy^/^ 



8nyUJ{X) 



-1 -2i 
-2i 1 



2'^rfiuj{\) \Ai 1 



-Ai 



+0 



nsm — 
2 



-3 



/(A)--3/2^-i(^) = / + Ai + A2 + O 



■ . 9^' 
nsm — 
2 



-3 



(108) 



A G dUs, 



where Ai and A2 denote the terms with nycJ(z) and n^u;(z), respectively. (Note that Ai(A) and 
A2(A) are analytic functions in Us \ {1} with poles of order 1 at A = 1.) This is an expansion in 
the inverse powers of nsin ^, and it holds uniformly for -^ < (/9 < vr, n > s, s > sq; and for A on 
the boundary dUs- 

Similarly, we define a conformal mapping for the neighborhood Us'- 

a;(A) = ln2(-<l>(A)), 

so that we have 



w(A) = -2usin2||l + ^ 



and 



cos 09 H — 

4 



u^O{i?)\, n = A + l, A G C/^, 



(109) 



^= -iV2u^/2sin-^(l + 0(u)), u = A + l. 

We have for the parametrix in Us'- 

P(A) = E[\)a^Q{n'uj[\))a^e-^^f^>'!{\r^''l\ A G Us. 



with 



E{\) 



1 



T2 



1 i 



N{X)f{Xr/' r M (7rn^^(A))'^3/2. 



Therefore, 
P(A)iV(A)-i = / + iV(A)/(A)'^-^/2 

+0 



-1 2i 
2i 1 



nsm — 
2. 



8ny^u;(A) 
/(A)--=^/2Ar-i(A) = / + Ai + A2 + O 



27n2tj(A) \-Ai 1 



1 4i 



nsm — 
2 



-3 



:iio) 



XedUs 



uniformly in A and in (p for — < ip < tt, n > s, s > sq. Similarly to the situation for Us, Ai(A), 
A2(A) in Us \ {—1} are analytic functions with poles of order 1 at A = —1. 
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5.3 Final transformation 

Let 



R{X) = S{X)N-\X), XgCMUsUUsU Si,2,3), 

R{X) = S{X)P-HX), XeUs\^i,2,3, (111) 

R{X) = S{X)P-\X), XeUs\^i 



Furthermore, set 



^1,2,3- 



R^V^-mV2. (112) 



It is easy to see that this function has jumps only on dUs, dUs, and parts of Si, and S3 lying 
outside of the neighborhoods Us, Us (we denote these parts S'j^ij*). Namely, 



i?+(A) = i?_(A)p-'^3^r(A) (^f^^^^il^^y2n J) N{X)-'V-^, X e S?|, 



(113) 



R+iX) = R^{X)V^^-^P{X)N{X)-'V^J, X G dUs 



R+iX) = i?_(A)P^'^3P(A)iV(A)-^PS, A G dUs, 

h. 

The jump matrices for R on dUs and dUs have the form 

/ + 0(p~^), p = nsin-, 

uniformly in A (as well as in a, (3 in compact sets) and in ip provided 

2s 

— < 99 < TT, n > s, s > sq. (114) 

n 

(A more detailed expansion is given by (108) and (110).) 

Let us now estimate the jump matrix on S^*^* for ip in the range (114) with sq sufficiently large. 
Below e' will stand for various positive constants independent of 93, n, and A. Denote x = KA, 
y = SJA, i.e., X = x + iy. Choose S°"* so that y is small, however, y > 5' sin ^ with some fixed 
6' > 0. Multiplying (83) with its complex conjugate and expanding in y gives: 

mx)f = i + ^JEi^4^^^ — AGsr. (115) 

V 1 — x^(sm ^ + X^ COS^ ^) 
Let |x| > e for some e > 0. Then we immediately obtain from (115): 

|c^(A)|2>l + e'sin|, 

for A G S°"*, and (p in the range (114). Now let |x| < e. We parametrise x as follows x = rsin^. 
Therefore < r < e/ sin ^ . If r > 1 , we have 

|<I>(A)|2 > 1 + -^ > 1 + - > 1 + e'sin^. 
1 + r^ r 2 
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On the other hand, if r < 1, we have, recalhng our condition y > 6' sin ^, 



Thus we conclude that the estimate |$(A)p > 1 + e'sin^ holds uniformly for A G S°"*, and if in 
the range (114). A similar estimate holds on S3"*. As N and / are bounded on A G S°H|' and, 
in particular, at A = 0, these estimates immediately imply that the jump matrix on S°'ij* can be 
written as 

/ + 0(e"^'''), p = nsin|, (116) 

uniformly for A G S^^g*, and 99 in the range (114). 

The above estimates for the jump matrices and the standard analysis of the R-RH problem 
(see [12]) imply that the R-RH problem is solvable for large sq a-i^d the solution has the form of the 
series: 

k-l 

R{X) = I + Y. RM) + 0(p-'), RM) = 0{p-'), (117) 

i=i 
uniformly for ip in the range (114), for q in a compact subset of the a-half-plane Jfta > —1/2, for 
/3 in a compact subset of the /3-plane, and for all A. 

Explicit expressions for Rk are obtained by collecting terms of the same order in the jump 
relations. Thus by (108) we have, in particular, that Ri, R2 satisfy the following RH problems. 
The functions i?i, R2 are analytic in C \ {dUs U dUs); 

i?i,+ (A) - i?i,_(A) = P^'^3A,(A)P^3^ 

_ _ _ _ (.lloj 

^2,+ (A) - i?2,-(A) = Ri,-{X)V^^-'Ai{X)V2 + V^^^A2{X)V2, XedUsU dUs; 

and i?i — 7- 0, i?2 — ^ as A — )■ 00. 

As Ai(A) is analytic in neighborhoods of the endpoints except for the simple poles at ±1, we 
have an expansion 

P^-3Ai(A)PS = 4^ + ^^'^ + ^^'^(A - 1) + 0((A - 1)2), as A ^ 1; 

V^--^A,{X)V2 = j^ + 5(2) + ij(3)(A + 1) + 0((A + 1)2), as A ^ -1, 

where A and B are constant matrices. It is now easy to verify that the RH problem for Ri has the 
following solution: 

I CT + m - -Di^AilAjCS. for A e Us U Us . 
For the function A2 we can write similarly: 



P^'^^^A2(A)PS = ^ + C^') + 0{X - 1), as A ^ 1; 
P^'^^A2(A)PS = ^ + ^^'^ + 0(A + 1), as A ^ -1, 



23 



:i2i) 



where C and D are constant matrices. A similar expression to (120) can now be written for i?2- 
However, we will need it only in the limit A — t- ±1. We then obtain: 

f \ , (122) 

i^2(-l) = -[A«, B«] + -(^(2)^(1) ^ ^(1)^(2) _ ^(D) ^ ^(2)2 ^ ^(3)^(1) _ ^(2)^ 



where [A, B] = AB - BA. We also obtain from (120), 



i2,(l) = ii^d) - A(2), i?i(-i) = -i^W - S(2). 



(123) 



5.4 Asymptotic form of the differential identity (49) 

We will now compute asymptotics of the r.h.s. of (49). First, we determine the components of the 
matrices A^^\ . . . , D^^> we need for the calculation below. First, the expansions for uj{X) in Section 
5.2 can be written as 



1 



1 



^/^) (2n)V2sin| 
1 1 



y^IJiX) -i(2n)i/2 



sinf 



(1 + LJi+U + UJ2+U^ + 0{U^)), 



{1 + OJi^U + UJ2-u'^ + 0{u^)), 



u = X — 1, 



u = X + 1, 



with 



^1 f 5 

^l± = ±2 I COSV9+- 



(The explicit expression for uj2± will not be needed below.) Therefore, expanding iV(A) we obtain 
from (108) 



1 



1 



16pV-* -1 

^°" 2 V-i(6 + 7o + 8Pi+) 



^(2) 



1 



70 



16p 
-i{6 + 70 - 8Vi 

-70 



^1+ I , ! ) + To 



-1 



(124) 



with 



p = n sm 



2' 



70 



4+4^?-- 



'!+• 



Here Pi-i- is given by (100). Furthermore, 



^(3) 



16p 



W2H 



1 



1 -i 
-i -1 

71 



+ wi+Fo + Pi 



-i(72 - 8P2 



with 



4 V -^(72 + 8P2+ + f ^f+) 
11 _o _ _ 4 



16t)3 ^ 
3 ^l+'l 



-71 



(125) 



71 = ^ - 2p2^ + 8P1+P2+ + -Vt+, 72 = 71 + 82??+ 



3. 



We obtain the matrices B^^', j = 1,2,3, by taking the matrices A^^', j = 1,2,3, and replacing p 
with —p; i with — i; wi,2+ with a;i^2-; ^1,2+ with X^i,2-; and the prefactor -^ of ^0 with —i^- 
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The components of C are also obtained from (108). For C (and similarly for D) we will only 
need the following combination of matrix elements: 



'-"21 ~ "-"12 



-3i 



Similarly, 



/I) 

"^21 



)(1) 
'l2 



2^2' 

-3i 

2V' 



^(2) _ ^(2) _ -3i 

"-"21 "-"12 - 25p2 



(2a;i++p2^). 



^(2) _ ^(2) 

^21 -^12 



3i 

2V 



(-2wi_ + Pil 



(126) 



(127) 



We now need to evaluate the asymptotics of trC+C_ and tr(7i(C-f -|-C_) from (49). Tracing 
back the transformations of the RH problem, we see that Y{z) for z close to Z4., so that A is close 
to 1 in the region ///, is given by 

Y{z) = 7"'^^roit;P^"'"3^ P(A) = ^(A)Q(n2w(A))e-"^/^^'^«/(A)~'"«/^ 

where E is given by (104), and Q by (107). Substituting this into (48) and expanding Bessel 
functions at ^ = 0, we obtain: 



c+ = ^^r^'nv^Rii) (J _\) Rii)-'v^--^T,w 



(128) 



Here we can write 

i?(i) = / + i?i(i) + ^2(1) + oip-^), ^(i)-i = / - ^i(l) - ^2(1) + Rjii) + oip-^). 

For definitiveness, we assume that the roots in (93) are chosen with the arguments from to 2tt, 
and that the point A corresponding to the limit in (48) approaches A = 1 from above. 

For C- we similarly have 

-i 1 



C. 



-2i 



r^'ToV'^Ri-l) 



1 



i?(-i)-ip^'^-^ro-S-«-3 



Therefore, we now easily obtain 

trC.C- 



1 + ti + t2 + O 



p- 



where (recall first (69), and then (123) and the above expressions for A, B) 

h = i(i?l,2l(l) - -Rl,12(l) - [^l,2l(-l) - i?l,12(-l)]) = -— (1 - Sin^ 

p \ 2 

The expression for t2 needs more work. First, we obtain it in the form 

t2 = i(i?2,2l(l) - i22,12(l) - [^2,2l(-l) " -R2,12(-l)]) + ^, 

where 



(129) 



(130) 



(131) 



(132) 



4J] = tr <Ri{l) 



Ri{l), 



+tr <^i2i(-l) 



iii(-l), 



-i 1 
1 i 

i 1 
1 - 



i 1 

1 -i 

-i 1 

1 i 



+ iT 



^i(l), 



i 1 



^i(-l), 



-i 1 
1 i 
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Denoting the matrix elements 

\ic —a J \-ic -a 

we obtain after a simple algebraic computation that 

17 = -(6 - c + 6 - c)^ + (a - a)^ - (c + c)(6 + b) - a{b + c) - a{b + c). 

Using again (123) and the above expressions for A^ B, we deduce from this formula that 

n = ^{A{Vl^ + Vl_) + 2(6^1+ - wi_) + 5 - 16(Pi+ + Vi^f). (133) 

On the other hand, we obtain from (122) and the above expressions for A, B, C, D: 

i(i?2,2l(l)-i?2,12(l)-[i?2,2l(-l)-i?2,12(-l)]) = ^(9- 22(l^i+ -L^i_) - 12(2??^ +Pt))- (134) 

Thus, substituting (133), (134) into (132) and then using the expressions for wi±, 'Df_^_, we obtain 
*2 = ^(1 - 6(wi+ - 6^1-) + 8Pi+Pi_) = ^ (a^ (l - sin |)' + /3' cos^ | - i cos^ ^\ . (135) 
Collecting together (131), (135), we finally have 

2 ■ 2 V5 „ • V? /-, . 'P\ 1 2'P , 2f. ■ ^\^ , «2 2 V? , ^/ -U 



tr (7+C_ = n sin 2an sin — ( 1 — sin — ) cos — \- a ( 1 — sin — ) + B cos \- 0(p 

+ 2 2V 2/4 2V 2/^ 2^'^ 

(136) 

We now turn our attention to the quantity trCi(C+ + C_). Recall that z = 1 corresponds to 
A = oo. We, therefore, immediately obtain 

Ci = aY{l)asY{l)-^ = a^^'-^T^a^T^^^'^'^K (137) 

In C+ + C_ we now need to take into account only the terms of order no less than Ri. Using (128), 
(129), and at the last step (131), we obtain 

trCi(C+ + C_) = an sin |(2 + ti) + 0(/9"^) = 2Qnsin|-2a2 ('i-sinl) +0{p'^). (138) 

Substituting this expression and (136) into (49), we finally have 
Proposition 9 ^s n — )■ cxd, 

d? ^ , , n2 , , 1-sin| i+4(a2_^2) ^/ 1 \ ^ ^ 

uniformly for ^ < (p < ir, n > s, s > sq, for some sq > 0, and uniformly in compact subsets of the 
half-plane 3?a > — 1/2 and of the plane /3 G C. 
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6 Proof of Theorem 1 

The asymptotic evaluation of the Toephtz determinant Dn{(p) is based on the integration of the 
differential identity (139) from cp to n — e with a small positive e. We have: 

{7r-e-^)ilnDny\n-e-lnDniTT-e)+lnDn{v)= f dot I{4>)d4>, (140) 



where I{(p) is the r.h.s. of (139). Fix n. Calculating the integral on the r.h.s. of (140), substituting 
for lnL'„(7r — e) the expansion (72) (and for (lnZ)„)'|7r-e its derivative), and taking the limit e — )■ 0, 
we obtain 



.2, ^ , o.„,^ , „,2m„ ^ , .■,^'P\ li„„ 1 



lni:'„(^) =nMncos- + 2(an + a^)lnf l + sin-j - -InnH In2 + 3C'(-1) 

/I ^ \ / 1 \ ^^ (141) 

-2a2ln2- ( - ^^ + a^]lnsm^ + O I ^-^] + (5„. 
\4 / 2 \nsm-|/ 

uniformly for — < ip < tt, n > s, s > sq, and where (5„ — )■ as n — )■ oo. This expansion is uniform 
in compact subsets of the a-half-plane Ka > —1/2 and of the /3-plane. Now substituting (38) and 
(141) into (30) and taking the limit n — t- oo, we obtain (8). D 

7 Bessel kernel. Proof of Theorem 4. 

In this section we set /? = 0. The weight (15) is then 

f{z,ip) = \z-l\'^"', z = e'\ (p<9<2TT-ip, (142) 

and f{z,(p) = on the rest of the unit circle. Note that / is an even function of the angle 6. Let 

oj{x,ip) = ^f. '^^ = 2"%^^, x = cos0. (143) 



sin 6*1 ^/l 



X 



This function is supported on [— l,cosc/?]. Consider the Hankel determinant with symbol uj{x,ip): 



COSip 



n— 1 



D^{ip) = det / x^+''uj{x,ip)dx 

^-1 / j,k=o (144) 

-I /"COS ip /"COS ip J^ 

= — / •••/ JJ (xj -Xfc)^ JJa;(xj,v9)(ia;j. 

There holds the following 

Lemma 10 Let K^essen ^^ ^^^ operator acting on L'^{0, (2s)^), with kernel (11). Then 



r,H(2s\ 

,(q-i/2)^ ^ " - lim " ^~' 

n^oo L>^(0) ■ 



det(/ - <eLi20L^(o,(2.)^) = jim -^^j^. (145) 
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Proof. First, as in the proof of the formula (33) of Lemma 6, one obtains that 

D^{ip) = D^{0) det(/ - Kn)mcos^,i) (146) 

where Kn is the operator acting on L^((cosc/?, l),dx) with kernel 

-X„__i Pn{x)Pn-l{y) - Pn~l{x)Pn{y) 



K^{x, y) = y/co{xMy) ^^""^ "^^^^ ^^^"^ "'^^-^^^ ''''' , (147) 

x„ x-y 

where Pk{x) = XkX^ + " " ") k = Nq, Nq + 1, . . ., with some No > 0, are the polynomials orthonormal 
on [—1,1] w.r.t. w{x,0): 

1 

Pk{x)x"^uj{x, 0)dx = x^^Skm, m = 0,l,... ,k. 
1 

The choice of the function (143) implies the following Szego relations between Pk{x) and the 
polynomials qk{z) = XkZ^ + " " "; given by (19) (see Lemma 2.5. of [10]): 

Pk{x) = ^ ^^^ ^ ^^{z-\2k{z) + z\2k{z'^)). (148) 

V27r(l + (?2fe(0)/x2fc) 

and for the leading coefficients, 



^ _^k^ . / i + g2fc(o)/x2fc 

Note that qn{z) = qn{z) as in our case /(e*^) = /(e~*^). 

Similarly to the proof of Lemma 6, we now set x = cos(2n/n), y = cos{2v/n), fix < u < s and 
< V < s and consider the limit of Kn{x, y) as n — )■ oo. We will now show that this double-scaling 
limit gives the kernel of K^^^^^^^ ■ 

First, it follows from Theorem 1.8. of [10] that 

.2„(0)_^A\ ^^1^0^' 



X2n \nj x„ 2 V^. 

Moreover, we have, 

X = cos (—] = 1 - ^ + o{n-^), z{x) = e^W", 
\ n J n^ 

and similarly for y = cos{2v/n). Therefore, by (143), 



Vio{x,OMy,0) = (-] {uvr-^/\l + 0{n-')). (149) 



28 



Moreover, taking zi = e^***'" G C+ and using the expression (21) of Lemma 5, we have 
q2{n+k){zi) = (2n)" ^^^^^^"U (l + a, 1 + 2a, 4iu{l + k/n)){l + 0{l/n)) 



(2nr^(l + «) 



r(l + 2a) 
fc = 0,-l. 



0(1 + a, 1 + 2a, 4m) -\ 4>'(1 + a, 1 + 2a, Aiu) 

n 



[l + 0{l/n)). 



Substituting this into (148), we then obtain an expression for P„(a:)P„_i(y) — Pn-i{x)Pn{y) in 
terms of the confluent hypergeometric functions and their derivatives at ib4m, ib4w. Removing 
the derivatives with the help of the standard relation 

(j)' {a, c, x) = —{4>{a + 1, c, x) — (pia, c, x)) 

X 

reducing then, by Kummer's transformation (27), the terms with the arguments —Aiu, —Aiv to 
functions of the arguments 4iu, Aiv, and making use of the following standard recurrence relation 

(c — a)(t){a — 1, c, x) + {2a — c + x)(j){a, c, x) — a(j){a + 1, c, x) = 0, 

we obtain 

Pnixi)Pn-lix2) " P„-l(xi)P„(x2) = " ^n^^o^l (^^)'"""''^"'"'^ " f ^ + ^ ( ')) 

7rr^(l+2a) n\ \nJJ 

X {{u + v)[(l){l + a,l + 2a,4m)(/)(a, 1 +2a,4w) - 0(1 + a,l + 2a,4m)0(a, 1 + 2a, 4m)] ^ ^ 
+{u- v) [0(1 + a, 1 + 2a, 4m)0(l + a, 1 + 2a, 4m) - 0(a, 1 + 2a, 4m)0(a, 1 + 2a, 4m)] } . 

Thus, in the difference Pn{x)Pn-i{y) — Pn-iix)Pniy) the main terms in n dropped out leaving the 
ones of order 1/n. Moreover, we did not need to know expressions for the terms 0(l/n) in Lemma 
5, as their contribution to the terms of order 1/n dropped out as well. 

We now employ the recurrence relations (36) to express 0(a, 1 + 2a, x) and 0(1 + a, 1 + 2a, x) 
in terms of 0(a,2a, x) and 0(1 + a,2 + 2a, x) and then use the connection with Bessel functions 
(4). We obtain recalling (149): 

~, ^ _ 'n? [uJa+l/2{2u)Ja_il2{2v) - vJa+l/2{2v)Ja~l/2{2u)] [^ , ^fl 

J^n\X,y) — 2 2 I '1 



which leads to 



1 



.(a-l/2)..„_x2 /o„,^2^ 



i^^ I -2^^ ) ^"(^'^) = -^Be7sL2{{2uY, {2vY), (151) 



where —Ksg^sssen ^^^^ "-"^ ((2-5)^)0) (note the reversed direction) and 



2{x-y) 



^BesseZ2(^' v) = o/^ A ' ' a = a - 1/2, (152) 



which, by the relation zJa+i{z) = aJa{z) — zJ'^{z), is equivalent to (11). 

The convergence of the determinants follows from the convergence of the kernels as in Lemma 
6, and we obtain the statement (145) from (146). D 
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We now evaluate the r.h.s. of (145). First, note that D„ (0) is a Hankel determinant whose 
symbol uj{x,0) (143) is supported on [—1,1] and has two Fisher Hartwig singularities at x = —1 
and X = 1. Therefore, the asymptotics of D^ (0) are given by a particular case of Theorem 1.20 
from [10]. Namely, 

G(l/2 + a) 

uniformly in compact subsets of the half-plane ?R-a > —1/2. 

In order to evaluate D^[2s/n), we use a connection formula between Hankel and Toeplitz 
determinants given by Theorem 2.6 in [10]. The formula is written in terms of the matrix elements 
of y(2n) (^^^) ^39) g^^^ fo^ ^ ^ 2s /n as foUows: 

The asymptotic expression for the Toeplitz determinant D2n {-^) is given (uniformly in com- 
pact subsets of the half-plane Ka > —1/2) by (141) with n replaced by 2n and with (p set to be 

2s/n: 

In D2n ( — ) = -2s^ + 4as - ( q^ + - j In s -F a^ In n - ( 2q2 + - 

1 /1\ (''') 

+-ln7r + 21nG(l/2)+0f-j +(5„(s), n ^ oo, 

where (5„(s) — )■ as n — )■ oo. 

It now remains to estimate Y-l^^^ (z) at z = —1,0,1. The A-images of these points are A = 
0, — ztan^,oo, respectively. All of them lie in the regions where Y^'^'^'[z) is approximated by the 
outside parametrix N{X). From the expressions (75,81,85-87,111,112,117), we obtain: 

y(^)(2n)^^2n.3^r(,tan|)-i(/ + PSO(/>-')^r=0^(M^))(5(^^) l)*W'"'^^ (156) 

which is valid for z in neighborhoods of z = —1,0, 1. In a neighborhood of —1 we assume that 
\z\ > 1 and then b{z) = /(z)~^^(z)~'^"z^"'. In neighborhoods of and 1, b{z) = 0. 

We will need the values of T>{X) at A = 0, — Jtan^,oo. Analyzing the integral in (94) with 
/3 = 0, we obtain 

Pl'aanl) =Z)('-itan|) = ('l + sinD", P(0) = 2". (157) 

From (95) we have 

/ 4sin^ \" 

2?oo = T— -^ • (158) 

V 1 + sm ^ / 



expression for the 11 element of (156): 



2 



Noting that in the definition of A'^(itan ^), a(itan ^) = e*''^ "^i/^^ we can write the following 



y/2«)(^) ^ i^2n^(^)2n f^^^^^^y v{\{z))-^ (e''^^-^'^/^a{\) + e''^^-^^ ' ^ a-\\) 



(159) 
+b{z)V{\{z))i-'^ (e*(^-")/^a(A) - e~^(^-")/^a-i(A))] (1 + 0{p-^)). 
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Since by (73), 



^(-1) = -1, *(0)=7-\ ^(l)=7-'(l + sm|), 



we obtain from (159) 

Yl^''\-l) = j^''2-''(l + sm^Y\cos^ + sm^] (l+Oip-^)) = 2-''[l+0{s-% n > s^ (160) 



2 



4 4 



where the second equation is obtained by substituting if = 2s/n, and considering n > s^ and s 
large. Similarly, 



and 



yr^(0)=sin^(l + O(p-i)) = -[l + O(s-i)], n>s, (161) 



/-. , ■ (£\2(n+a) 2s a. 

^^^a) = ^^S^^°^^(^ + ^(^")) = 2if^tl + 0(.-)], n>s\ (162) 

Note that (160-162) are uniform in compact subsets of the half-plane Kq > —1/2. 

Substituting (160-162) and (155) into (154) and taking the square root, we obtain 

D^ n±\ = vr"+3n^2-(--i)'+i"-"'+iG(l/2)s-H"-i)%-^'+(2"-i)«[l + 0{s-% n > s\ 

(163) 
The branch of the square root is fixed by the fact that D2n > for a G M, and by the uniformity 
of the asymptotic expansion in a. Finally, substituting (163) and (153) into (145), we finish the 
proof of Theorem 4. D 

8 Appendix 

Here we show that the operator i^("'^) with kernel (1) on L^(— s,s), where Ka > —1/2 and 

1/2 

a lb /? 7^ —1, 2, . . ., is trace class. Note first that since gj (x) is bounded on {—s, s), it is sufficient 
to show that the operator K with the following kernel 

K{x, y) = 1^1 1^1 [e'^y-^'^n + a + /3, 1 + 2a, 2ix)Ml + a - /3, 1 + 2a, -2iy) 

x-y (164) 

_(.i(^-y)^{l j^ a + 13,1 + 2a, 2i7/)0(l + a - /3, 1 + 2a, -2ix)] 
is trace class. Expanding the confluent hyper geometric functions in series, 

oo oo 

e"*''(^(l + a + /3, 1 + 2a, 2ix) = ^ /U^a;", e*''(/)(l + a - /3, 1 + 2a, -2ix) = ^ A^x", 

n=0 m=0 

where fin, An are determined using (2), we can write the kernel (164) in the form: 



xUym _ yU^m 



k{x, y) = Ixn^r E E l^n\ra— "-^ ■ (165) 

^-^ ^-^ X — y 

n=Om=0 
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We will now show that the trace norm in Lp'{—s, s) 

x-y 
for some C > 0. Together with the straightforward estimates 

t 



n 

A'n ) ^n^ — T) ^ ^ 1, 

n! 



for some 6 G M, the inequality (166) implies that K, and hence K^'^'^' , is trace class. 
To prove (166), set first m = 0, n > 1. Then we have for some (7 > 



„n „,n n 1 



IIfI — z — 1^1 lli<2^llly| IIl2(-s,s)IIfI IIl2(-s,s) - 

^ l^2(fe + Ka) + lJ (^2(n-fc + KQ)-lJ " ^ ^ 

n— 1 



(^y^g|{2(A:+3f?a)+l)+|(2(n-A:+gf?a)-l) ^ (^^gn+23?a 



A;=0 
which gives (166) for m = 0, n > 1. If m > 1, n > 1, we can assume that n > m and write 

ri.n„.m „.nrf,m r^n—m „.n~m 

I l^x (/ y X I iQ, I ,ci+m'^ y I \a+m 

y — y 5 

X — y a^ — y 

which is of the same form as (167) with a replaced by a + m, and n, by n — m. Hence, we complete 
the proof of (166). 
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